Delay times and reflection in chaotic cavities with absorption 
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Absorption yields an additional exponential decay in open quantum systems which can be described by shift- 
ing the (scattering) energy E along the imaginary axis, E + ih/2r a . Using the random matrix approach, we 
calculate analytically the distribution of proper delay times (eigenvalues of the time-delay matrix) in chaotic 
systems with broken time-reversal symmetry that is valid for an arbitrary number of generally nonequivalent 
channels and an arbitrary absorption rate t^ 1 ■ The relation between the average delay time and the "norm- 
leakage" decay function is found. Fluctuations above the average at large values of delay times are strongly 
suppressed by absorption. The relation of the time-delay matrix to the reflection matrix S'S is established at 
arbitrary absorption that gives us the distribution of reflection eigenvalues. The particular case of single-channel 
scattering is explicitly considered in detail. 
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There is a growing interest in statistical properties of the 
Wigner-Smith matrix Q{E)=-i%S^dS/dE with S(E) 

being the scattering matrix at the collision energy E, in the 
cases of chaotic scattering and transport in disordered media 
131 . In the resonance scattering, the matrix element Q cc r de- 
scribes the overlap of the internal parts of the scattering wave 
functions in the incident channels c and c' This di- 

rectly relates the Wigner-Smith matrix to the effective non- 



Hermitian Hamiltonian Tt = H 



'iVV' of the unstable in- 



termediate system as follows (henceforth 7i=l) | 



Q(E) = V 



1 



1 



[E-Uy E — Ti 



V . 



(1) 



The Hermitian part H stands here for the closed counterpart of 
the system while the amplitudes V£ describe the coupling be- 
tween N interior and M channel states. The random-matrix 
theory approach is usually adopted to simulate the compli- 
cated intrinsic motion |5JeIlZI|- 

The known analytical results H H E El E El Q are 
restricted to the idealization neglecting absorption. The latter 
is, however, always present to some extent under laboratory 
conditions, being one of the sources of a coherence loss in 
quantum transport. This has dramatic consequences for the 
statistical observables flallal . Necessity of proper account- 
ing of finite decoherence 1 1711 was recently emphasized jT3l in 
order to remove a certain discrepancy between theory I19t 12011 
and experiment [16] on conductance distributions in quantum 
dots. Reflection in a weakly absorbing medium turned out to 
be directly related fl5l, EH l22ll to the time-delay matrix with- 
out absorption. Recent experiments 12311 in microwave cav- 
ities demonstrated that the absorption (due to the skin effect 
in the walls) may be strong, leading to an exponential decay 

In this paper we show that representation Q in terms of the 
effective Hamiltonian allows us to extent the consideration to 
the case of an arbitrary absorption. The nature of the exponen- 
tial decay caused by absorption can be easily understood from 
the following model consideration which actually goes back 
to the concept of the spreading width in nuclear physics 12411 ; 
see 1 4] for the recent developments. In addition to coupling 



to continuum (scattering) states the originally closed system 
is considered to be also coupled to the background compound 
environment. The latter has a very dense spectrum with the 
mean level spacing Ab g being much smaller than the corre- 
sponding one A of the closed system, Ab g <C A. When 
the coupling strength v 2 is large enough to mix background 
states, v 2 > A 2 ^, the original levels acquire the damping or 



spreading width Tj = 2itv /Ab g B4l 12411 . Corrections to the 
resulting exponential decay show up at the time t* ~ A,^ 1 
and, therefore, can be safely ignored on mesoscopic scale of 
the Heisenberg time tn = 2ir/A <C t* we are interested in. 
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FIG. 1: An open cavity with absorption in walls modelled by cou- 
pling to (a) the compound background or (b) an infinite number of 
fictitious channels with vanishing transmission probabilities. 

In the absorption limit of continuous spectrum of the back- 
ground, when an irreversible decay into walls takes place, this 
description becomes equivalent to that achieved in the frame- 
work of the Biittiker's model of dephasing in mesoscopic con- 
ductors; see Fig. 1 l25ll . One considers 1171 fl^l fictitious 
scattering channels in addition to M real ones. The vanish- 
ing transmission T$ — * of the fictitious channels is assumed 
to be compensated by their large number — ► 00, the di- 
mensionless absorption rate 7 = M^T^ being kept fixed 1 20] . 
Then the anti-Hermitian part of the effective Hamiltonian H, 
which describes coupling to (all) open channels, splits readily 
as X^reai V n V m* + S n m ?a HU into the escape contribution 
(first term) and damping one with T a = t^ 1 = 7A/271". An 
associated with the last term exponential decay lasts up to the 
characteristic time t* = t#/ y/jT^ 1 26] being large as com- 
pared to tn ■ 

The consideration presented suggests that nonzero absorp- 
tion is equivalent to the purely imaginary shift E + |T a = E 1 



2 



of the energy in the Green's function (E — H)^ 1 of the open 
system without absorption as long as resonance scattering 
far from the channel thresholds is concerned ll27ll : see also 
1I22I Eijll . This is in agreement with available data on correla- 
tions of S matrix elements in cavities with absorption 1 23 ] . 

In what follows we consider the time-delay matrix with ab- 
sorption Q 7 = Q(Ery), with Q from treating -f — T a tH as 
a phenomenological parameter. The important relation for the 
reflection matrix 



r — SySy — 1 r a Q-j 



(2) 



follows directly from the definition of the scattering matrix 
S 7 = S{E 7 ) = 1 - iV^{E 1 - UY X V, which is subunitary 
(R < 1) at nonzero absorption. This relation gives Q 1 the 
meaning of the matrix of unitarity deficit and generalizes lim- 
iting expressions of Refs. I2lll22ll valid at weak absorption to 
the case of arbitrary T a . Q 1 isaMxM Hermitian, positive- 
definite matrix and, therefore, has real positive eigenvalues q c , 
the so-called proper delay times. They were recently studied 
in much detail for the case of zero absorption II 3lll4l . Even 
a weak absorption modifies their statistical properties signifi- 
cantly, as will be shown below. 

We begin with the calculation of the average total delay 
time g t ot = Qi + • • • + Qm = tr Q 7 , where the bar denotes 
the ensemble average. Making use of the invariance of the 
trace under cyclic permutations and the following relation 
VVt = i[(E 7 - W)t - {E 1 - U)\ - r o , one gets 



gtot = Im Tr 



e^-u 



-T a Tr 



1 



1 



E y - H (E 1 - W)t 



(3) 



where Tr acts in the TV-dimensional intrinsic space of res- 
onances. The first term is known I& |9J to be equal to the 
Heisenberg time tn ■ To calculate the second one, it is instruc- 
tive to go to the time domain and to exploit the well-known 
relation between the Green's function and the time evolution 
operator cxp(-ifit). This enables us to represent l|3} in the 
following form: 



Ttot 



gtot 



= i-r„ 



dtt 



(4) 



where P(t) = {l/N)TT(e iUH e- im ) is the "norm-leakage" 
decay function introduced in Ref. 1 26] . The average delay 
time within the cavity becomes smaller due to additional dis- 
solution in the walls. The average weighted-mean reflec- 
tion coefficient (r) = M~ l trR is correspondingly given by 
(r) = 1 - 77t t/M. 

P{t) can be calculated by means of Efetov's supersymme- 
try technique |5, 29], which becomes now a standard analyt- 
ical tool. Here we only state the corresponding result for the 
case of preserved time-reversal symmetry (TRS): 



P(t) 



1 p OO pOO 

dX / dXi / dX 2 
-1 Ji Ji 

M r 



ti*i)s(- 



:/(Ai)II 



t A1A2-A 
*H 2 

{9c + A) 2 



(.9c + A 1 A 2 )2-(A 2 -l)(A2-l) 



where fi(Xi) = (1 - A 2 )/(A 2 + A 2 + A 2 - 2AAiA 2 - l) 2 
and f(Xi) = (2A 2 A^ - X\ - X 2 2 - A 2 + l)/4. The quantities 
g c = 2/T c — 1 > 1 are related to the transmission coeffi- 
cients T c = 1 — I Sec 1 2 JH, which determine the openness 
strength of the system (without absorption), referring T = 1 
(0) to the completely open (closed) one. For reader's conve- 
nience, we note that the result for the case of broken TRS |26] 
follows from by removing there the A2 integration and set- 
ting A2 = 1 everywhere in the integrand save the integration 
measure /i(Aj) = (Ai — A)~ 2 in this case 1 30 ] . It is worth also 
pointing out the relation between P(t), Eq. (|5}, and the au- 
tocorrelation function of the photodissociation cross section 
l3lll . The exact (in the RMT limit N — ► 00) equation @ is 
valid for any symmetry and will be also derived below using 
a different way. 

The "norm-leakage" is identical to unity when the system 
is closed (hence the norm). Its time dependence is solely due 
to the openness of the system and has been thoroughly studied 
in 1 26] that allows us to understand the qualitative dependence 
of cjftot on absorption. The typical behavior P(t) ~ IIc=i ^ + 
(2//3)T c i/iff] _,3/2 , with 0=1(2) standing for preserved (bro- 
ken) TRS, is the simple exponential exp(— t J2 C Tc/tn ) at 
small enough times. In the so-called "diagonal approxima- 
tion" j2(J, which neglects the nonorthogonality of the res- 
onance wave functions and becomes asymptotically exact at 
large t, P(t) turns out to be related by the Laplace transform 
P diag (i) = f o °°dr e- Ft p(T) = (e~ rt ) r to the distribution 
p(T) of resonance widths. One gets readily from @ that 
Ttot = (r/(r + r a )) r within this very approximation. The 
simple interpolation formula r to t ~ (1 + 7/ Ylc ^c) 1 with 
corrections of the order of min[l/7, l/Se^e] becomes ex- 
act as the absorption rate 7 and/or the total (dimensionless) 
escape width ^ c T c grows. 

We proceed further with an analysis of the distribution of 
the proper delay times V{q) — M^ 1 J2 C ^(l ~ 1c)- For the 
sake of simplicity, we restrict ourselves to the case of broken 
TRS (the unitary symmetry class). The factorized representa- 
tion ([Hi of Q 1 enables us to use the same method developed in 
Ref. 1 14] to treat the zero absorption case. Thus, we skip all 
standard technical details, indicating only essential ones. As 
usual, the jump of the resolvent G(z) = M _1 tr (z — Q 7 ) _1 
on the discontinuity line along q = Kcz > determines the 
seeking distribution as follows: "P(q) = 7r _1 Im G(q — iO). 
Due to the factorized structure of Q 1 , G(z) can be then rep- 
resented in the form suitable for subsequent supersymmetry 
calculation J32il . We find the following expression for the de- 
termining part K{C = z/t H ) = M( 2 (t H G(z) - 1/C): 



K(0 = 1 



dX 1 



dX 



M 

n 



9c + A 



( d 



, ,MAi)A(A) 



(6) 



) Here b^X^e^W^^^Ioif^/il - 7CHA 2 - 1)] and 
/ 7 (A)=e-( 1 -7C/2)-A/C Jo[ ,y(i _ 7C)(1 _ A2)]> with h{x) 

(5) [>^o(^)] being the modified (usual) Bessel function. The re- 
solvent G(C) given by Eq. is an analytical function of 



3 



the complex variable £ for the negative values of Re £ and, 
therefore, can be expanded there in Taylor's series. One 
finds directly from the definition of G that tnG(z) = 1/C + 
tr Q 7 / (M^Hh) + ■ ■ ■ for large (, relating thus q w to the co- 
efficient of the second term of this expansion. On the other 
hand, this coefficient is given just by K(— oo) which can eas- 
ily be calculated from l|6} to reproduce exactly equation @. 

An analytical continuation in (|6j to the region of positive 
t = Re £ requires more care as compared to the case IU4I1 
of zero absorption, where it was achieved by a proper defor- 
mation of an original integration contour. First we make the 
following decomposition in partial fractions: 



1 



M 

n 



9c + A 



1 



M 

E 



i 



Ai - A g c + Ai Ai - A g a + M 



n 



9b + A 



The contribution from the term (Ai — A) -1 leads to an exact 
cancellation with the first term in (|6). [This is not surprising 
since the product term in 0, the channel factor, which deter- 
mines solely the strength of system openness, reduces at Ai=A 
to unity, resulting G(z) — l/z in this case identically.] The in- 
tegration over Ai gets completely decoupled from that over A 
in the contribution from the rest sum. Making use of the table 
integrals |33J, one finds that 



D d\ ie 



9 + A] 



-I (a 



A? 



1) 



dp- 



-m>, 



(7) 

with shorthands s = t _1 — 7/2 and a = t — T7. Just 
this term @ has a nonzero imaginary part, thus the distribu- 
tion, at positive t — iO. A close inspection of the r.h.s. of 
Eq. @ shows, that the imaginary part is determined by the 
integration region s — a < p < s + a, resulting at the end in 
irIo(ay/ g 2 — l)0(r _1 — 7), with the step function Q{x). 

We arrive finally at the following general expression for the 
probability distribution of the proper delay times: 



1 



M 



M ^ 

c=l 
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- Pi b ^ 

OVx) 



for < r < 7 1 , and V{t) = otherwise. Here 

1 



B,. 



loiviaVgT 



9a - 9c 



(8) 



(9a) 



+1 



dX 



e~ vsX J (zWT 



A 2 ) J] (g a + A) . (9b) 



The obtained result is valid for arbitrary absorption strength 
and arbitrary transmission coefficients of M generally 
nonequivalent channels. The limit of zero absorption 11411 
is correctly reproduced. The case of statistically equivalent 
channels can easily be worked out by performing the limiting 
transition g c — > g = 2 /T — 1 for all c. At last, the distribution 
function P R (r) = M^Y^J^ 
ues r c = 1 



[r — r c ) of reflection eigenval- 
11c/ tH follows readily from (|8j as 



P R (r) =j- 1 V[j- 1 (l~r)] , 0<r<l. 



(10) 



We see that the absorption rate 7 enters the distribution in 
a highly nontrivial way. This is expected to be true for any 
distribution function and is contr asted with a correlation func- 
tion of, say, S matrix elements S* ab (E 1 )S a >b'{E 1 + e). The 
corresponding form-factor 1 23] (the Fourier transform of the 
correlation function) differs from that 01 of the zero absorp- 
tion case simply by the presence of an additional exponential 
term e -7 *' tff . The most striking effect of finite absorption on 
the time-delay distribution is likely to consist in suppression 
of the universal long-time tails T -M0/2-2 (T^, EH [T^, [3] 
at t > 7 _1 1 34]. To understand this fact qualitatively, we note 
that the delay time q(E) « T n /[(E-E n ) 2 + \(T n + T a ) 2 } in 
a vicinity of a given resonance with the energy E n . The max- 
imal value of this single-resonance contribution is attained at 
E = E n , being q max = AT n /(T n + T a ) 2 < 1/T a for any 
value of the (positive) escape width T n . 

We analyze now the important case of single-channel scat- 
tering, AI = 1, in more detail. The explicit expression to be 
obtained from (|8jl reads as follows: 



V(r) = 



h{ot\fi 



7 2 

r2 \ v .'/- - l)(cosh - - - sinh ^) 



2 N 

2' 



-- sinh - [gsl (a^g 2 - 1) 
7 2 

-a^g^lhia^/g^l)}) 



(11) 



We have explicitly checked the normalization of this distribu- 
tion to unity and verified relation @ for the first moment. This 
function should be compared to the more simple expression 
7> 7=0 (t) = T-^d/d^e-s^Ioir-^g 2 - 1) valid at 
zero absorption. Figure 2 shows the behavior of "P(t) in two 
limiting cases of the weakly and perfectly open system. One 
sees in the first case that the maximum of the distribution func- 
tion at the small time r ~ (2j) _1 « T/4< 1 gets more pro- 
nounced and narrow as the absorption rate 7 grows. At larger 
values of r the distribution is exponentially suppressed with 
P(7~ 1 ) r; (7 2 /T)e~ 7/T . The latter is contrasted with the 
behavior in the case of perfect coupling, g = T = 1, when 
Vt=i{t) = T- 2 e- 1 / T [l + (1 - r)(e 7 - l)/r 7 ] and V^- 1 ) 
could be rather large. We relate these distinctions to pecu- 
liarities in fluctuations of the resonance escape widths in the 
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FIG. 2: (color online). The distribution Jl H of the time delay in 
single-channel scattering for different values of the absorption rate 7 
at weak (T = 0.1, left) and perfect (T = 1, right) coupling. 
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FIG. 3: (color online). The reflection coefficient distribution in the 
single-channel cavity at four experimental realizations [36] of the 
absorption rate and transmission coefficient (see the text for details). 
The values (27, T) correspond to I: (0.56, 0.12), II: (2.42, 0.75), III: 
(8.4, 0.98), and IV: (48, 0.99). 

two cases considered. The width distribution p(T), which is 
known exactly 1 10] for any T and M, has the simple expo- 
nential form e~ ViH l T when coupling is small, T < 1, and 
the power law behavior ~ T~ 2 at T > i^ 1 when T = 1. The 
ratio of the widths V ~ T a determining V(t ~ 7 _1 ) is, there- 
fore, exponentially small in the first case and only a power in 
the second. This conclusion holds for any finite M. 

The sharp border at r = 7 _1 of the obtained distribution 
is the direct consequence of equation (0 with the absorption 
rate fixed to a constant. Although, as shown above, the value 
7 3 (7~ 1 ) of the jump may be exponentially small when cou- 
pling is weak, a generic exponential suppression should be 
intuitively expected at large values of delay times r 3> 7 _1 - 
Indeed, for the time St a wave-packet oscillating in the cavity 
with the frequency A/2ir on average experiences (A/27r)<5i 
collisions with the walls, yielding the probability Tj,(A/2iv)8t 
to be absorbed into one of fictitious channels. The total 
reflection is then estimated as R ~ (1 - T0(A/27r)<5i) M *, 
giving e~ 7<5 */ tff in the limit of fixed 7 = M^T^ as — > 00 
and — > 0. It is instructive, therefore, to define alterna- 
tively through the following relation R = e~~ Ta ® R the ma- 
trix Qn, which we call the matrix of reflection time-delays. 
The positive definite matrix Qr is related to Q 1 as Qr = 
— r^ 1 ln(l — V a Qi) that leads to the following connection 

V R (T r ) = e-^V[y- 1 (l-e-^)], r r > , (12) 

between the corresponding distributions 7 3 fl(r r ) and V{f) of 
proper delay times (eigenvalues of Qr and Q 7 , respectively). 
Both Q r and Q 1 reduce to the same Wigner-Smith matrix Q 



in the limit of vanishing absorption. The difference between 
them becomes noticeable at finite 7. Still both distributions 
coincide up to the time appreciably less than 7 _1 . They start 
to differ at larger times, when V{t) has the cutoff whereas 
'Pr{ t 2> 7 _1 ) e~ 7T is exponentially suppressed. 

Finally, we discuss the distribution Pr(t) of the reflection 
coefficient r = |5 7 | 2 = 1 — 7T in the single-channel cavity. 
This distribution at arbitrary values of 7 and T is explicitly 
given by Eqs. dlOl and il Q , reproducing exactly the recent 
result [35] obtained by a different method. In the particular 
case of perfect coupling it simplifies further to the expression 
P T=1 (r) = (l-r)- 3 e- 7 /( 1 - r )[7(eT-l) + (l+7-eT)(l-r)] 
found earlier I22ll . For the case of preserved TRS ((3 — 1), the 
reflection coefficient distribution in a microwave cavity has 
recently been measured 113611 . Our distribution Pij(r) at the 
values of absorption and transmission realized in this experi- 
ment under compulsory (although not surprising in the RMT) 
rescaling 7 to 7/3/2 with /? = 1 is shown on Fig. 3. (This cor- 
responds to replacing our parameter 7 in (II 01 and il Q with 
T w /2 of Ref. B36I1 .) That should roughly take into account the 
difference between the symmetry class of our analytical result 
(f3 = 2) and that of the experiment. Such a replacement is 
expected to become more efficient as absorption grows. The 
trend is clearly seen from the distribution sharply peaked near 
7* ~ 1 at weak absorption (7 <C 1) to the Rayleigh distribu- 
tion P R (r) ~ (7/3/2)e- r7 ' 3 / 2 |28], see also |22], reproduced 
correctly at strong absorption (7 ^> 1) and perfect coupling 
(T = 1). Figure 3 is in good qualitative agreement with 
the experimental data reported in Ref. 1 36] (see Figs. 4 and 
6 there), which becomes even quantitative as absorption gets 
stronger. The rigorous analytical treatment for the case of pre- 
served TRS is still lacking, being under current investigation. 

In summary, we have calculated the general distribution 
of proper delay times and reflection coefficients in an open 
chaotic system (e.g., billiard) with broken TRS at arbitrary ab- 
sorption. Finite absorption leads to strong suppression of fluc- 
tuations at large values of delay times, making the distribution 
narrower around the mean. The latter as well as the mean re- 
flection coefficient are found to be related to the "norm leak- 
age" decay function. The particular case of single-channel 
scattering is paid appreciable attention, when discussion of 
available experimental data is also given. 
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No. 03-02-16151 (D.V.S.) is acknowledged with thanks. 



[1] E. Wigner, Phys. Rev. 98, 145 (1955). 
[2] F. T. Smith, Phys. Rev. 118, 349 (1960). 

[3] For a general review of the concepts of the time delay and its 
applications see C. A. A. de Carvalho and H. M. Nussenzveig, 
Phys. Rep. 364, 83 (2002) as well as Ref. C3]. 

[4] V. V. Sokolov and V. Zelevinsky, Phys. Rev. C 56, 31 1 (1997). 



[5] J. J. M. Verbaarschot, H. A. Weidenmiiller, and M. R. Zim- 

bauer, Phys. Rep. 129, 367 (1985). 
[6] C. W. J. Beenakker, Rev. Mod. Phys. 69, 731 (1997). 
[7] Y. Alhassid, Rev. Mod. Phys. 72, 895 (2000). 
[8] V. L. Lyuboshitz, Phys. Lett. B 72, 41 (1977); Yad. Fiz. 27, 948 

(1978) [Sov. J. Nucl. Phys. 27, 502 (1978)]. 



5 



[9] N. Lehmann, D. V. Savin, V. V. Sokolov, and H.-J. Sommers, 

PhysicaD86, 572 (1995). 
[10] Y. V. Fyodorov and H.-J. Sommers, J. Math. Phys. 38, 1918 

(1997); Phys. Rev. Lett. 76, 4709 (1996). 
[11] V. A. Gopar, P. A. Mello, and M. Buttiker, Phys. Rev. Lett. 77, 

3005 (1996). 

[12] Y. V. Fyodorov, D. V. Savin, and H.-J. Sommers, Phys. Rev. 

E 55, R4857 (1997); D. V. Savin, Y. V. Fyodorov, and H.-J. 

Sommers, Phys. Rev. E 63, 035202(R) (2001). 
[13] P. W. Brouwer, K. M. Frahm, and C. W. J. Beenakker, Phys. 

Rev. Lett. 78, 4737 (1997); Waves Random Media 9, 91 (1999). 
[14] H.-J. Sommers, D. V. Savin, and V. V. Sokolov, Phys. Rev. Lett. 

87, 094101 (2001). 
[15] E. Doron, U. Smilansky, and A. Frenkel, Phys. Rev. Lett. 65, 

3072 (1990). 

[16] A. G. Huibers, S. R. Patel, C. M. Marcus, P. W. Brouwer, C. I. 

Duruoz, and J. S. Harris, Jr., Phys. Rev. Lett. 81, 1917 (1998). 
[17] M. Buttiker, Phys. Rev. B 33, 3020 (1986). 
[18] E. R. P. Alves and C. H. Lewenkopf, Phys. Rev. Lett. 88, 

256805 (2002). 

[19] H. U. Baranger and P. A. Mello, Phys. Rev. B 51, 4703 (1995). 
[20] P. W. Brouwer and C. W. J. Beenakker, Phys. Rev. B 55, 4695 

(1997) [Erratum: Phys. Rev. B 66, 209901(E) (2002)]. 
[21] S. A. Ramakrishna and N. Kumar, Phys. Rev. B 61, 3163 

(2000). 

[22] C. W. J. Beenakker and P. W. Brouwer, Physica E 9, 463 (2001). 
[23] R. Schafer, T. Gorin, T. H. Seligman, and H.-J. Stockmann, J. 

Phys. A: Math. Gen. 36, 3289 (2003). 
[24] A. Bohr and B. R. Mottelson, Nuclear Structure (Benjamin, 



New York, 1969), Vol. 1. 

[25] It is worth noting that the case of absorption, when dissolution 
in walls occurs, should be contrasted with that of dephasing, 
when the particle number in the cavity is conserved; see 1 20]. 

[26] D. V. Savin and V. V. Sokolov, Phys. Rev. E 56, R491 1 (1997). 

[27] Only in this case the smooth energy dependence of the coupling 
amplitudes can be neglected, so that any scattering quantity de- 
pends explicitly on E by means of (E-H)" 1 only. 

[28] E. Kogan, P. A. Mello, and H. Liqun, Phys. Rev. E 61, R17 
(2000). 

[29] K. B. Efetov, Supersymmetry in Disorder and Chaos (Cam- 
bridge University Press, Cambridge, UK, 1996). 

[30] This is the usual procedure in the crossover regime of gradually 
broken TRS; see 1 12, 29]. Note also a misprint in 1 26]: /(Ai) = 
(Ai — A 2 )/4 is correct. 

[31] Y. V. Fyodorov and Y. Alhassid, Phys. Rev. A 58, R3375 
(1998). 

[32] The generating function (4) of 1 14] can be again written in terms 
of the 27V x 2N determinants Det(A±), where now A± = 

- {zg 1 - §r o ) - ^Vl -z±r n Vi - i{E-H)a 3 . 

[33] j™dte- st I l±aVt?=I] = e^ s2 - a2 /Vs 2 - a 2 for s > 
a > 0. For the calculation of the imaginary part we use: 

J^dt e~ a9t cos(aVT~^¥) /VT~ 7 ~t 2 ~ = tt/o^Vs 2 ~ -1L 
[34] Similar suppression was previously observed in Ref. |21] for 

the time-delay distribution in single-channel reflection from 

one-dimensional amplifying random media. 

[35] Y. V. Fyodorov, cond-mat/0304673 

[36] R. A. Mendez-Sanchez et al. cond-mat/0305090 



